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, We apply the general results of the kinetic theory of systems with long-range interactions to 

particular systems of physical interest. We consider repulsive and attractive power-law potentials of 
' interaction r ' with 7 < d in a space of dimension d. For 7 > 7c = (d — l)/2, strong collisions must 

, be taken into account and the evolution of the system is governed by the Boltzmann equation or 

» I ■ by a modified Landau equation; for 7 < 7c, strong collisions are negligible and the evolution of the 

' system is governed by the Lenard-Balescu equation. In the marginal case 7 = 7c, we can use the 

Landau equation (with appropriately justified cut-offs) as a relevant approximation of the Boltzmann 
and Lenard-Balescu equations. The divergence at small scales that appears in the original Landau 
Xf^ ' equation is regularized by the effect of strong collisions. In the case of repulsive interactions with a 

neutralizing background (e.g. plasmas), the divergence at large scales that appears in the original 
Landau equation is regularized by collective effects accounting for Debye shielding. In the case of 
attractive interactions (e.g. gravity), it is regularized by the spatial inhomogeneity of the system and 
its finite extent. We provide explicit analytical expressions of the diffusion and friction coefficients, 
and of the relaxation time, depending on the value of the exponent 7 and on the dimension of space 
d. We treat in a unified framework the case of Coulombian plasmas and stellar systems in various 
dimensions of space, and the case of the attractive and repulsive HMF models. 
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I. INTRODUCTION 



The dynamics and thermodynamics of systems with long-range interactions is currently a topic of active research 
in physics [11-01 • A system is said to be long-ranged if the potential of interaction decays at large distances as 
P5 , r^'' with 7 < d. In that case, the potential energy u{r)r'^~^ dr diverges at large scales implying that all the 
O ' particles interact with each other and that the system displays a collective behavior. For such systems, the mean 
field approximation becomes exact in a proper thermodynamic limit N — >■ -l-oo In 0, Q (Papers I and II), we 
have developed a general framework to tackle the kinetic theory of systems with long-range interactions and we have 
presented the basic kinetic equations (Vlasov, Landau, Lenard-Balescu, Fokker-Planck). These equations were initially 
I> • introduced in plasma physics [§-[l2j and stellar dynamics [l3l - [l5j but they are actually valid for a much larger class 
'tJ ■ of systems with long-range interactions. In this paper, we consider specific applications of this general formalism. 
Wc provide explicit analytical expressions of the diffusion and friction coefficients, and of the relaxation time, for 
plasmas and stellar systems in various dimensions of space, and for the attractive and repulsive Hamiltonian Mean 
• , Field (HMF) models. Our approach provides a unified framework to treat these different systems. In this paper, we 
■ restrict ourselves to systems that are spatially homogeneous or for which a local approximation may be implemented. 
The case of spatially inhomogeneous systems is more complicated and must be treated with angle-action variables as 

in 

We first consider repulsive and attractive power-law potentials of interaction r~'' with 7 < d in a space of dimension 
d, generalizing the traditional Coulombian and Newtonian potentials (corresponding to 7 = d — 2). For A ^ 1, where 
A is the effective number of particles^, the evolution of the system is dominated by weak collisions and collective 
^ , effects. We can therefore make a weak coupling approximation and expand the equations of the BBGKY hierarchy 
in powers of 1/A (see Appendix IXj) . For A — > +00, we get the Vlasov equation. At the order 1/A, the evolution of 
the distribution function is governed by the Lenard-Balescu equation. However, for singular potentials, the Lenard- 
Balescu equation may present divergences at small scales reflecting the importance of strong collisions that have been 
neglected in the weak coupling approximation. For power-law potentials, there exist a critical index 7c = (d — l)/2 
(for Coulombian or Newtonian interactions, it corresponds to a critical dimension dc = 3). For 7 > 7c (i.e. d > 3 
for Coulombian or Newtonian potentials) , strong collisions must be taken into account while collective effects may be 
neglected. In that case, the evolution of the system is described by the Boltzmann equation or by a modified Landau 



^ For attractive interactions, A = n\j represents tlie number of particles in the Jeans sphere, which is of the order of the total number of 
particles N. For repulsive interactions with a neutralizing background, A = n\jj represents the number of particles in the Debye sphere. 
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equation (see Appendix |B]) . The relaxation time^ scales as ~ A^'^~'^~^^/'^tD, where tu is the dynamical time. For 
7 < 7c (i.e. d < 3 for Coulombian or Newtonian potentials), strong collisions are negligible while collective effects are 
important. In that case, the Lenard-Balescu equation is rigorously valid. The relaxation time scales as ^ Md. 
For 7 = 7c (i.e. d = 3 for Coulombian or Newtonian potentials), strong collisions and collective effects must be taken 
into account but their importance is weak. The Boltzmann equation is marginally valid provided that a large-scale 
cut-off is introduced at the Debye length. Similarly, the Lenard-Balescu equation is marginally valid provided that a 
small-scale cut-off is introduced at the Landau length. In that marginal situation, we can use the Landau equation 
(with appropriately justified cut-offs) as a good approximation of the Boltzmann and Lenard-Balescu equations. The 
relaxation time scales as tn ~ (A/lnA)i£). These scalings agree with those obtained by Gabrielli et al. [2l| based on 
an extension of the Chandrasekhar binary collision theory. The binary collision theory is usually adapted to systems 
with short-range interactions. It can be used when 7 > 7c and it is marginally valid when 7 = 7c (corresponding 
to 3D Coulombian or 3D Newtonian potentials). For 7 < 7c it produces a divergence at large scales. By contrast, 
the Lenard-Balescu theory is adapted to systems with long-range interactions. It can be used when 7 < 7c and it is 
marginally valid when 7 = 7c. For 7 > 7c it produces a divergence at small scales. Therefore, these two theories are 
complementary to each other. 

We discuss in detail the physical rcgularization of the divergences that occur in the Landau equation. The small- 
scale divergence is regularized by the effect of strong collisions. They are taken into account in the Boltzmann 
equation or in the modified Landau equation. For repulsive potentials with a neutralizing background (as in plasma 
physics), the large-scale divergence is regularized by collective effects accounting for Debye shielding. In a plasma, 
since the Coulomb force between like-sign (resp. opposite-sign) charges is repulsive (resp. attractive), the sign of 
the polarization cloud surrounding a test charge is opposite to that of the test charge. As a result, the interaction is 
screened on a distance of the order of the Debye length A^i [l^]. Collective effects are accounted for in the Lenard- 
Balescu equation. This equation docs not present any divergence at large scales, contrary to the Landau equation, 
and the Debye length appears naturally. For attractive potentials (as in stellar dynamics), the large-scale divergence 
is regularized by the spatial inhomogeneity of the system and its finite extent. The Jeans length Aj [2^, which is 
of the order of the system's size, represents a natural large-scale cut-off. Spatial inhomogeneity is accounted for in 
the Landau and Lenard-Balescu equations written in angle-action variables (l6l - l20| . These equations do not present 
any divergence at large scales, contrary to the case where we make a local approximation. In the gravitational case, 
collective effects lead to a form of anti-shielding. The test star draws neighboring stars into its vicinity and these add 
their gravitational force to that of the test star itself. The "bare" gravitational force of the test star is thus augmented 
rather than shielded. The polarization acts to increase the effective gravitational mass of a test star. As a result, 
collective effects tend to increase the value of the diffusion coefficient and reduce the relaxation time of the system 

[Ullllll. 

For 3D plasmas, the Lenard-Balescu equation is marginally valid provided that a small-scale cut-off is introduced 
at the Landau length in order to take into account the effect of strong collisions. The relaxation time scales as 
(A/lnA)tD where A is the number of electrons in the Debye sphere. The Landau equation with a small-scale cut-off 
at the Landau length and a large-scale cut-off at the Debye length provides an excellent approximation of the Lenard- 
Balescu equation. For 2D plasmas, the Lenard-Balescu equation is rigorously valid. The relaxation time scales as 
At£). The Landau equation with a cut-off at the Debye length is a poor approximation of the Lenard-Balescu equation 
but the discrepancy is not dramatic. In particular, the velocity dependence of the diffusion coefficient is the same, 
up to a numerical factor. Therefore, we may use the Landau equation by "adapting" the value of the large scale 
cut-off (it is a fraction of the Debye length). For ID plasmas, the Lenard-Balescu equation is rigorously valid but it 
trivially reduces to zero. The relaxation time of the system as a whole scales as h'^to while the relaxation time of 
a test particle in a thermal bath scales as Md. The Landau equation leads to wrong results. It predicts that the 
diffusion coefficient of a test particle in a thermal bath decays like a Gaussian while in reality, when collective effects 
are properly accounted for, it decays as \v\~^ . 

For 3D stellar systems, the Lenard-Balescu equation and the Landau equation written with angle-action variables 
are marginally valid provided that a small-scale cut-off is introduced at the Landau length in order to take into 
account the effect of strong collisions. The relaxation time scales as {N/\nN)to where N is the number of stars in 
the cluster. The Landau equation based on a local approximation with a small-scale cut-off at the Landau length and 
a large-scale cut-off at the Jeans length is relatively accurate. For 2D stellar systems, the Lenard-Balescu equation 
and the Landau equation written with angle- action variables are rigorously valid. The relaxation time scales as Nto- 



^ Here, we consider the relaxation time of a test particle in a thermal bath. For d > 1, it coincides with the relaxation time of the system 
as a whole. However, for spatially homogeneous systems in d = 1, the collision term vanishes and the relaxation time of the system as 
a whole is larger than Ktjj (see discussion in Paper H). 
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The Landau equation based on a local approximation produces a linear divergence at large scales but it may be used 
provided that a large scale cut-off is introduced and properly "adapted" (it is a fraction of the Jeans length). For 
ID stellar systems, the Lenard-Balescu equation and the Landau equation written with angle-action variables are 
rigorously valid. The relaxation time scales as Ntu- The Landau equation based on a local approximation leads to 
wrong results. In particular, it predicts a relaxation time of the system as a whole scaling as N'^to while the right 
scaling is Nto when spatial inhomogeneity is accounted for. 

For the HMF model, explicit expressions of the diffusion and friction coefficients, and of the relaxation time, can 
be obtained in the homogeneous phase. The relaxation time of the system as a whole is larger than NtD (scaling 
presumably as N'^to) while the relaxation time of a test particle in a bath scales as Ntn- The diffusion coefficient 
depends on the temperature. For the attractive HMF model, there is a critical point {Ec,Tc) = (3/4,1/2) and the 
diffusion coefficient increases close to the critical point implying a decrease of the relaxation time. 

Some results of kinetic theory have been obtained previously in particular cases, using different approaches, and 
we give the corresponding references. However, the originality of our approach is to develop a general framework to 
treat the kinetic theory of systems with long-range interactions (Papers I and H) and, from this framework, consider 
particular systems. This provides a more unified description of systems with long-range interactions. 

The paper is organized as follows. In Sec. [Hi we consider power-law potentials of interactions and introduce 
appropriate notations. In Sec. IIIIl we provide an estimate of the relaxation time for systems interacting with power- 
law potentials and we discuss how it depends on the value of the exponent 7. In Sec. IIVI and fVl we treat the case of 
Coulombian plasmas and self-gravitating systems in d dimensions. In Sec. IVIl we consider the attractive and repulsive 
HMF models. The Appendices gather additional results of kinetic theory that complete the main discussion of the 
paper. 

II. POWER-LAW POTENTIALS OF INTERACTION 

In this section, we consider power-law potentials of interaction of the form r"'' with 7 < d in a space of dimension d. 
We introduce appropriate notations that generalize those commonly used in the study of plasmas and self-gravitating 
systems. 

A. Repulsive interactions: plasmas 

The potential of interaction of a one-component Coulombian plasma with a neutralizing background in d-dimensions 
is the solution of the Poisson equation Au = — 5'rf(e^/77i^)(5(r) where Sd is the surface of a unit sphere, — e is the charge 
of the electron, and m is its mass. Its Fourier transform is (27r)''{t(fc) = SdC^ / {'m?k'^)- 

More generally, we consider a repulsive power-law potential whose Fourier transform may be written as (27r)'*u(fc) = 
Sde^ /m?k'^~°' . The potential of interaction in physical space is of the form u{y) = Kd.aiSde^ / rri?) , where 
Kd,a = J dqe'i Vg^-" = Tr^'i/'^2°'-'^r{{d -2 + a)/2)/r((2 - a)/2) is a constant depending on d and a (its 
precise expression is not needed since only the Fourier transform of the potential of interaction matters in the kinetic 
theory). The exponent of power-law decay is 7 = c? — 2 -I- a. The Coulombian potential corresponds to a = 0, i.e. 
7 = d — 2. For a = 2 — d (i.e. 7 = 0), u(r) = Kd{e'^ /m?) In |r| is the logarithmic potential. 

We introduce the parameter (3 = l/{mv'^) where Vm is a typical velocity of the system (for a system at statistical 
equilibrium, 13 = l/ksT is the inverse temperature). We define the Debye wavenumber ko and the plasma pulsation 
wp by the relations = Sde^Pp/m and Wp = Sdpe^kf^/inn? where p = nm is the mass density. They are related to 
each other by ojp ~ ko / \/ 13m — kuVm. We also introduce the dynamical time to ~ i-iJp^ = \/]3rn/kD = ^/{koVm) = 
^ I \/ Sdpe^k'^ j m? . It represents the typical time needed by a particle to travel the Debye length A_d ~ k'^ with 
the velocity Vm ^ {/3m)~^/'^. Finally, we introduce the parameter A = nXf) that gives the number of electrons 
in the Debye sphere. The normalized potential 'i'i{k) = —{2ttY u{k)l3mp introduced in Paper II may be written as 
r?(fc) = -(fci5/fc)2-". 

The number of charges in the Debye sphere scales like A ~ „(2-a-d)/(2-a)(y/g2)d/(2-a)^ p^^. ^y^q Coulombian 
potential (a = 0), we get A ^ 7i(2-<^)/2y<i/2^gd^ ^j^g number of charges in the Debye sphere is always an increasing 
function of the temperature. By contrast, its dependence with the density depends on the dimension of space. In 
d = 3, A ^ n~^/^T^/^/e^ decreases with the density. In d = 2, A ^ T/e^ does not depend on the density. In d = 1, 
A ^ n^/^T^/^/e increases with the density. 

Remark 1: in the previous formulae, e represent the charge of the electron only when a = 0. However, when a 7^ 0, 
it is always possible to write the potential of interaction in the above form (where e may be regarded as a coupling 
constant) so as to facilitate the connection with Coulombian plasmas when a = 0. 
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B. Attractive interactions: self-gravitating systems 

The potential of interaction of a self-gravitating system in d-diniensions is the solution of the Poisson equation 
Au = SiiG6{r) where G is the gravitational constant. Its Fourier transform is (27r)'^M(fc) = —SdC/k"^. 

More generally, we consider an attractive power-law potential whose Fourier transform may be written as 
{2'KYu[k) = —SdG/k'^~". The potential of interaction in physical space is u(r) = Kd.aSdG /\r\'^'~^'^" . The expo- 
nent of power-law decay is 7 = d—2 + a. The Newtonian potential corresponds to a = 0, i.e. 7 = d — 2. For a = 2 — d 
(i.e. 7 = 0), u{r) = KdG In |r| is the logarithmic potential. 

We introduce the parameter /3 = l/{mv'^) where Vm is a typical velocity of the system (for a system at statistical 
equilibrium, /3 = l/ksT is the inverse temperature). We define the Jeans wavenumber kj and the gravitational 
pulsation ujg by the relations fcj~" = SdG/3mp and ujq = SdGpkj where p = nm is the mass density. They are 
related to each other by ujg = kj / j3m = kjv„i. We also introduce the dynamical time to = = \fW^I^J ^ 
l/{kjVm) = ^/ yJSdGpk". It represents the typical time needed by a particle to travel the Jeans length Aj ^ fcj^ 
with the velocity Vm ~ (^m)~^/^. Finally, we introduce the parameter A = nXj that gives the number of particles in 
the Jeans sphere. The normalized potential r](k) = —{2TrYu{k)Pmp may be written as r]{k) = {kj / k)'^^" . 

For attractive power-law potentials in d dimensions with a 2 — d, the virial theorem reads 2K + (d — 2 + a)W = 
where K is the kinetic energy and W the potential energy (see Appendix I of (26l|). The kinetic energy may be 
estimated by K ^ Nmv'^^ and the potential energy by \W\ ^ N^rn^G / R'^^^'^" where R is the system's size. Since the 
virial theorem implies K ~ \W\, we find from the foregoing expressions that /3GMm/ R'^~'^~^°' ~ 1. When a = 2 — d, 
the virial theorem reads 2K — GM^J2 = leading to the exact result (w^) = GM/2 (the velocity dispersion can take 
a unique value in a steady state) [2a]. At statistical equilibrium K = dNksT /2, so we get (3GMm ~ 2d. Using 
p ~ M/R"^, we find in all cases that Aj ^ R. Therefore, the Jeans length Aj represents the typical size R of the 
system. Accordingly, the parameter A — nX'j ^ N gives the total number of particles in the system. 

Remark 2: in the previous formulae, G represent the gravitational constant only when a = 0. However, when a ^ 0, 
it is always possible to write the potential of interaction in the above form (where G may be regarded as a coupling 
constant) so as to facilitate the connection with self-gravitating systems when a ~ 0. The statistical mechanics of 



systems with attractive power-law interactions has been studied in 27 1 



C. Characteristic lengtlis 

There are three characteristic lengths in the problem: 

• The length / ^ n~^/'^, where n is the numerical density, gives the typical distance between particles. 

• The Debye length Ad = (m/S'^e^/Jp)^/^^""-' for repulsive potentials gives the effective range of the interaction 
due to the screening by opposite charges, and the Jeans length Aj ~ (l/S'^G/Jmp)-'^/^^""' for attractive potentials 
gives the typical size R of the system. These expressions rely on a mean field approximation so they are valid only 
for long-range interactions (7 < d i.e. a < 2). We shall denote commonly the Debye and Jeans lengths by A^ and 
define A = nAf . 

• The Landau length Al is the distance at which binary collisions become "strong" . This corresponds to an impact 
parameter yielding a deflexion at 90°. This happens when the energy of interaction e^/A''^^+" or G'm^/A'^~^+" 
between two particles becomes comparable to their kinetic energy mv'^^. Equating these two quantities, we get 
Al ~ (e^/mv^^)^/'''~^+") and Al ^ {Gm/v^y^^'^~^~^°'\ These definitions make sense only for decaying potentials 
(7 > i.e. a > d — 2). Introducing the notations defined in the preceding sections, the Landau wavenumber may be 

rewritten as k^ ^ (n/fc^^")"'^^'"'^ ~ A-'^/^'^^^+^^fc^. For attractive interactions, using Xj ^ R and A ~ iV, we get 

Al ~ i?/iVi/(''-2+a). 

Combining the previous expressions, we obtain 



I J V Al / V Al 



A. (1) 



For 2 - d < a < 2 (i.e. < 7 < d) and A > 1, we find that Al < ^ < A^. 

The coupling parameter F = Epot/Euin is equal to the ratio between the potential energy Epot ^ e^/n"^'^"^"'"")/'^ or 
Epot ~ Grr? I n~'^'^~'^^'^^ I of two particles separated by the average inter-particle distance I ~ 71^^/'' and the typical 
kinetic energy of a particle E^^in ^ mv^ = 1//3. This yields F = e'^n^'^~^~^°'^/'^ /mv^ or F = Gmn'^'^"'^^"-^/'^ /v"^. Using 
the previous relations, we find that F ^ l/A'^^""^/''. The weak coupling approximation corresponds to F ^ 1. For 
long-range interactions (7 < d i.e. a < 2), the conditions A ^ 1 and F ^ 1 are equivalent. We may also define the 
couphng parameter g = E*^^/Ekin where E*^^ ^ e^/A^^^^" or E*^^ ^ Gm'^ / Xj~^~^" is the potential energy of two 
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particles separated by the Debye length or by the Jeans length. In that case, we find g ~ 1/A. The parameter g is 
called the plasma parameter or the graininess parameter. It is equal to the reciprocal of the number of particles in 
the Debye or Jeans sphere. 



III. ESTIMATE OF RELAXATION TIME FOR SYSTEMS WITH POWER-LAW POTENTIALS 

The normalized Fourier transform of the power-law potentials defined previously may be written as 



^7^ -^[-f) ' (2) 

where e = +1 for attractive interactions and e = —1 for repulsive interactions (in the repulsive case, we assume 
that there exist a neutralizing background that maintains the spatial homogeneity of the system). The Coulombian 
potential (plasmas) corresponds toe=— 1,q; = and ks — kjj (Debye wavenumbcr) and the Newtonian potential 
(gravity) corresponds to e = +1, a = and fcs = kj (Jeans wavenumbcr). 

To estimate the relaxation time, we consider the Landau equation (11-5)^ in which strong collisions and collective 
effects are neglected (see Appendix |X] for more details about the approximations made in this section) . Therefore, we 
only consider weak collisions. We also assume that the system is spatially homogeneous or that a local approximation 
may be implemented. In the thermal bath approach, the diffusion coefficient is given by Eq. (11-45) with 

The relaxation time may be estimated by t^*'' ~ v'^^/'Dd (see Paper II). Defining the dynamical time by ~ ui'^^ ~ 
(vmks)^^ and performing the change of variables k = k/kg, we obtain 

'n' - ,+00^,. to, (4) 

where A = nk~'^. For repulsive interactions A — nk^'^ represents the number of particles in the Debye sphere and 
for attractive interactions A = nfcj'' ~ N represents the total number of particles in the system. The relaxation time 
scales like Ktu except if the integral in Eq. ^ diverges. In that case, it must be regularized at small or large scales 
and the regularized integral may depend on A. We must distinguish different cases according to whether a (or 7) is 
larger, smaller, or equal to the critical index 

S-d / d-l 



(5) 

For a Coulombian or a Newtonian potential (a = 0,7 = d— 2), this critical index corresponds to the critical dimension 
dc = 3. 



A. The case ac < ct < 2 



If ttc < a < 2 (i.e. 7c<7<rford>3 for a Coulombian or a Newtonian potential), the integral ([3]) converges 
for fc — )• (large scales) implying that collective effects (and spatial inhomogeneity effects for attractive interactions) 
are weak. On the other hand, it diverges algebraically for k +00 (small scales). This divergence is regularized by 
taking strong collisions into account. Heuristically, we can introduce a small-scale cut-off at the Landau length at 
which collisions become strong. Performing the integrals in Eqs. ([3]) and (U) with k^ax kj^ = A^i-s+a kg, we obtain 

Vd = ^ — r ^1^, tR"" - KT^to. (6) 

2(27r)'^-3 d-3 + 2a R d W 



^ Here and in the following, (Il-n) refers to Eq. (n) of Paper II. 
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We now develop a more precise description (see Appendix IXj) . In the limit A ^ 1, the evolution of the system is 
dominated by weak collisions and collective effects (and spatial inhomogeneity for attractive interactions) . Three-body 
collisions are negligible. In principle, strong collisions should also be negligible but the divergence of the diffusion 
coefficient shows that they are important at small scales. Since the divergence is strong (algebraic), the Lenard- 
Balescu equation is not applicable. On the other hand, collective effects and spatial inhomogeneity do not seem to 
be crucial since the diffusion coefficient ([3]) converges at large scales. If we only take weak and strong collisions into 
account, we obtain the Boltzmann equation. This equation does not present any divergence. Since weak collisions 
dominate over strong collisions when A ^ I, we can expand the Boltzmann equation for small deflexions (or directly 
use the Fokker-Planck equation). In that case, we obtain the modified Landau equation (jBI[) of Appendix [BJ This 
equation does not diverge at small scales since strong collisions have been accounted for. We note, however, that a 
more rigorous kinetic equation should take into account collective effects (and spatial inhomogeneity for attractive 
interactions) even if their presence is not required to make the integral converge at large scales. 



B. The case a < Oc 



li a < ac (i.e. 7 < 7c or d < 3 for a Coulombian or a Newtonian potential), the integral ([3]) converges for k — > +00 
(small scales) implying that strong collisions are negligible. On the other hand, it diverges algebraically for fc — > 
(large scales). In the case of repulsive interactions (e.g. plasmas), the divergence is regularized by collective effects 
which account for Debye shielding. Heuristically, we can introduce a cut-off at the Debye length and take fc,„i„ ^ kjj. 
In the case of attractive interactions (e.g. gravity), the divergence is regularized by the finite extent of the system. 
Heuristically, we can introduce a cut-off at the Jeans length and take kmin ~ kj. Performing the integrals in Eqs. 
and (21) with kmin kg, we obtain 

Ti — ^ ^ '"m^s bath a , 

We now develop a more precise description (see Appendix For repulsive interactions (e.g. plasmas), in the limit 
A ^ 1 the system is dominated by weak collisions and collective effects (strong collisions and three-body collisions can 
be neglected). Therefore, its evolution is rigorously described by the Lenard-Balescu equation (II-3). This equation 
does not present any divergence. In the thermal bath approximation, the diffusion tensor is given by Eq. (11-39) with 

We note that the large-scale divergence that occurs in the diffusion coefficient ([3]) is regularized by Debye shielding. 
If we make the Debye-Hiickel approximation, or consider small velocities |v| — > 0, the diffusion tensor is given by Eq. 
(11-45) with^ 

T,DH 1 ''^m^D f^°° k a 



If we neglect collective effects and introduce a large-scale cut-off at the Debye length, we get the Landau equation 
(II-5). In the thermal bath approach, the diffusion tensor is given by Eq. (11-45) with Eq. ([7]). However, the 
Landau equation with a cut-off at the Debye length is not quantitatively correct since the results depend strongly 
(algebraically) on the precise value of the large-scale cut-off. It may, however, provide a reasonable approximation of 
the Lenard-Balescu equation provided that the cut-off is suitably adapted to the situation. 

For attractive interactions (e.g. gravity), in the limit A f the system is dominated by weak collisions, spatial 
inhomogeneity, and collective effects (strong collisions and three-body collisions can be neglected). Therefore, its 
evolution is rigorously described by the Lenard-Balescu equation written with angle-action variables (l6l - [20| . If we 
neglect collective effects, it reduces to the Landau equation written with angle-action variables. These equations do 



The value of the integral is 
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not display any divergence. The large-scale divergence that occurs in the diffusion coefficient ([3]) is regularized by the 
finite extent of the system. If we neglect collective effects, make a local approximation, and introduce a large-scale 
cut-off at the Jeans length, we get the Vlasov-Landau equation (11-16). In the thermal bath approach, the diffusion 
tensor is given by Eq. (11-45) with Eq. ([7|). However, the Vlasov-Landau equation with a cut-off at the Jeans length is 
not quantitatively correct since the results depend strongly (algebraically) on the precise value of the large-scale cut- 
off. It may, however, provide a reasonable approximation of the Landau equation written with angle-action variables 
provided that the cut-off is suitably adapted to the situation. Finally, we note that collective effects tend to reduce 
the relaxation time (I7]-b) as explained in Appendix [Cj 



C. The marginal case a = Uc 

If a = ac (i-e. 7 = 7c or d = 3 for a Coulombian or a Newtonian potential), the integral ([3]) diverges logarithmically 
at small and large scales. This implies that both strong collisions and collective effects (or spatial inliomogeneity 
for attractive interactions) must be taken into account. However, their influence is weak since the divergence is only 
logarithmic. We are therefore in a marginal situation. Heuristically, we can introduce a small-scale cut-off at the 
Landau length'^ and a large-scale cut-off at the Debye length (for repulsive interactions) or at the Jeans length (for 

attractive interactions). Performing the integrals in Eqs. ([3]) and (|4]) with k^ax ^ = A"^^ kg and k„iin ^ kg, we 
obtain 

-n 1 2 vik. In A .^^^ A 

We now develop a more precise description (see Appendix IX)) . For repulsive interactions (e.g. plasmas), in the limit 
A ^ 1 the evolution of the system is dominated by weak collisions and collective effects. Three-body collisions are 
negligible. In principle, strong collisions should also be negligible but the divergence of the diffusion coefficient ([3]) 
shows that they are important at small scales. Therefore, the evolution of the system is rigorously described by Eqs. 
(jAlI a) and (jAlI b) with T = 0. These equations do not present any divergence. However, it is difficult to make them 
more explicit. The usual strategy is to consider successively the contribution of collisions with small and large impact 
parameters, and then connect these two limits. If we ignore collective effects and introduce a large-scale cut-off at 
the Debye length Xmax ~ ^d, we get the Boltzmann equation. This equation does not diverge at small scales since 
strong collisions arc taken into account. This equation is marginally valid since the divergence at large scales is weak 
(logarithmic). For A 3> 1, weak collisions dominate over strong collisions and we can expand the Boltzmann equation 
for small deflexions (or directly use the Fokker-Planck equation). In the dominant approximation InA S> 1, we get 
the Landau equation (II-5) with the diffusion coefficient (jlip with In A = lii{Xmax / ^l) in which the Landau length Xl 
appears naturally (see Appendix [D|. If we ignore strong collisions and introduce a small-scale cut-off at the Landau 
length Xmin = Xl, we get the Lenard-Balescu equation. This equation is marginally valid since the divergence at 
small scales is weak (logarithmic). In the thermal bath approach, the diffusion tensor is given by Eq. (11-39) with 



We note that the large-scale divergence that occurs in the Landau diffusion coefficient ([3]) is regularized by Debye 
shielding. If we make the Debyc-Hiickcl approximation, or consider small velocities |v| — > 0, the diffusion coefficient 
is given by Eq. (11-45) with^ 



10 [K(l + d)/2 + l]2 



^ We assume d > 1 because the Landau length is ill-defined in d = 1 when a = cic 
^ The value of the integral is 



L 



In 1 + Ad-i 



\d-l 

d + 1 

1 + 



(13) 



with kL/kn = A2/('^-i). 
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In the dominant approximation In A ^ 1, the Lenard-Balescu diffusion coefficient T>d{x) and the Debye-Hiickel 
diffusion coefficient T)^^ can be approximated by Eq. ([TT]) with InA = ln(A£i/Amm) in which the Debye length 
appears naturally. In conclusion, in the dominant approximation In A ^ 1 , the evolution of the system is rigorously 
described by the Landau equation with a small scale cut-off at the Landau length and a large-scale cut-off at the 
Debye length. These cut-offs are not ad hoc but they are justified by the above arguments. 

For attractive interactions (e.g. gravity), in the limit A ^ 1 the evolution of the system is dominated by weak 
collisions, collective effects, and spatial inhomogeneity. Three-body collisions are negligible. In principle, strong 
collisions should also be negligible but the divergence of the diffusion coefficient (jS]) shows that they are important at 
small scales. Therefore, the evolution of the system is rigorously described by Eqs. (|A2l -a) and (|A2l -b) with T = 0. 
These equations do not present any divergence. However, it is difficult to make them more explicit. Again, the 
usual strategy is to consider successively the contribution of collisions with small and large impact parameters, and 
then connect these two limits. If we ignore collective effects and introduce a large-scale cut-off at the Jeans length 
Amaa; = Aj, we get the Boltzmann equation. This equation does not diverge at small scales since strong collisions 
are taken into account. This equation is marginally valid since the divergence at large scales is weak (logarithmic). 
For A ^ 1, weak collisions dominate over strong collisions and we can expand the Boltzmann equation for small 
deflexions (or directly use the Fokker-Planck equation). In the dominant approximation In A 3> 1, we get the Vlasov- 
Landau equation (11-16) with the diffusion coefficient (|lip with InA = ln(A„,ax/AL) in which the Landau length A^ 
appears naturally (see Appendix [D|. If we ignore strong collisions and introduce a small-scale cut-off at the Landau 
length Xmin = Al, we get the Lenard-Balescu equation written with angle-action variables [i6, .20j. This equation is 
marginally valid since the divergence at small scales is weak (logarithmic). If we neglect collective effects, it reduces 
to the Landau equation written with angle-action variables. These equations do not diverge at large scales since they 
take into account the finite extent of the system. If we neglect collective effects and make a local approximation 
(which is justified by the fact that the divergence of the diffusion coefficient at large scales is weak), we get the 
Vlasov-Landau equation (11-16). This equation presents a logarithmic divergence at large scales (fc — >■ 0). Since this 
divergence is cured by the finite extent of the system, it is natural to introduce a large-scale cut-off at the Jeans scale. 
Since this procedure is essentially heuristic, the precise value of the large-scale cut-off is not known. However, in the 
dominant approximation In A 3> 1, the precise value of the numerical factor is not important since it appears in a 
logarithmic factor and we finally obtain the Vlasov-Landau equation (11-16) with the diffusion coefficient pT|) with 
In A = ln(A,//A„iin)- In conclusion, in the dominant approximation In A ^ 1, the evolution of the system is reasonably 
well described by the Vlasov-Landau equation (11-16) with a small scale cut-off at the Landau length and a large-scale 
cut-off at the Jeans length. This equation is not exact (the rigorous equation is the Lenard-Balescu equation written 
in angle-action variables with a small-scale cut-off at the Landau length justified by the previous arguments) but it 
may be used as a simplified kinetic equation. Finally, we note that collective effects tend to reduce the relaxation 
time (jlll b) as explained in Appendix [Cl 

D. Existence of quasi stationary states 

The previous scalings of the relaxation time agree with those obtained by Gabrielli et al. [2l| from the Chandrasekhar 
binary collision theory^. For a < ac, the Lenard-Balescu equation is rigorously valid and the relaxation time scales as 
Ai£). For a = ac, the Lenard-Balescu equation is marginally valid and the relaxation time scales as (A/ In A)i£). This 
corresponds to the usual situation encountered in plasma physics and stellar dynamics (a = 0) since the dimension of 
space d = 3 is critical for Coulombian and Newtonian interactions. For a > ac, the Lenard-Balescu equation is not 
valid and the relaxation time scales as A<'-2+°t£). 

For a < 1 (corresponding to 7 < d — 1), the relaxation time diverges when A — > +00 (for a < ac, the relaxation 
time increases linearly with A and, for ac < a < 1, it increases less rapidly than A). In that case, the Vlasov 
regime becomes infinite in the thermodynamic limit and quasi stationary states (QSS) may form as a result of violent 
relaxation [l^. On the other hand, for a > 1 (corresponding to 7 > d — 1), the relaxation time tends to zero when 
A — >■ -|-oo. In that case, as emphasized by Gabrielli et al. |2l| . no QSS can form (except if the potential includes a 
sufficiently large soft core). In the intermediate case a = 1, the relaxation time is independent on A. 

Remark 3: If we define the mean free path by A '--^ VmtR, we find that A - \(l~o^)/(d~2+a)^^ fQj. 

ac < a < 2, 

A ^ (A/ In A)As for a = ac, and A ~ AA^ for a < ac- For a < 1. wc find that X ^ \s when A 1 so that the mean 



Their approach assumes that collective effects and spatial inhomogeneity can be neglected. Our approach is more general since it can 
take these effects into account as explained previously. 
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free path is much larger than the Debye length or the Jeans length. This is another manner to understand why the 
Vlasov equation becomes exact when A — >■ +00. 

Remark 4-' For specific power-law potentials, the dynamics may be particular. For 7 = —2 (harmonic pot ential) 
there is no (violent or slow) relaxation and the system oscillates permanently as discovered by Newton |29l - [31| . 



IV. COULOMBIAN PLASMAS 



For Coulombian plasmas, the Fourier transform of the potential of interaction is 



i^-y^k)-^. m--^, (15) 



and the Fourier transform of the Debye-Hiickel potential is 



Q,p2 1 1.2 
(2-)W(fc) = ^^^-^, VDHik)^-j^^. (16) 



In physical space, we have u{r) = {e'^/m?)r ^ and uuni^) = (e^/m^)e '"'o''/r in c? = 3, u{r) = — (e^/rn^) ln(r) and 
UDH{r) = (e^/TO^)A'o(fc£)r) in d = 2, u{x) = — (e-^/m^)|a;| and udh{x) — (e^/m^fc£))e~*''°l^l in d = 1. 



A. 3D Coulombian plasmas 



If we neglect strong collisions and collective effects, the evolution of the system is described by the Landau equation 
(II-5). In the thermal bath approach, the difi'usion tensor is given by Eq. (11-45) with 

2(27r)i n k ^ ^ 

The integral diverges logarithmically at small and large scales. This implies that both strong collisions and collective 
effects (Debye shielding) must be taken into account. However, their effect is weak since the divergence is logarithmic 
(we are in the marginal case of Sec. IIII Cp . If we introduce a small-scale cut-off kmax = = Afci) at the Landau 
length at which binary collisions become strong and a large-scale cut-off fc„„;„ = kjj at the Debye length at which the 
interaction is shielded, we get 

^ 1 InA , , 

1^3 = — , (18) 

where A = nfc^"^ 3> 1 represents the number of electrons in the Debye sphere. In plasma physics, InA is called the 
Coulomb logarithm. 

Strong collisions are taken into account in the Boltzmann equation. This corresponds to the two-body encounters (or 
impact) theory. In this equation, the integral over the impact parameter does not diverge at small scales. Therefore, 
the correct treatment of strong collisions regularizes the logarithmic divergence that appears at small scales in the 
Landau equation. However, for a Coulombian potential in d = 3, the integral over the impact parameter diverges 
logarithmically at large scales. Therefore, the Boltzmann equation is marginally valid provided that a large-scale 
cut-off is introduced at the Debye length Xmax ~ (see the next argument). If we expand the Boltzmann equation 
for small deflexions (or directly use the Fokker-Planck equation), and consider the dominant approximation In A ^ 1, 
we obtain the Landau equation (II-5) with the diffusion coefficient (|18l) with In A = lii{Xmax / ^l) in which the Landau 
length Al appears naturally (see Appendix |D|) . 

Collective effects are taken into account in the Lenard-Balescu equation (II-3). This corresponds to the wave theory. 
In this equation, the integral over the wavenumber does not diverge at large scales. Therefore, the correct treatment of 
collective effects regularizes the logarithmic divergence that appears at large scales in the Landau equation. However, 
for a Coulombian potential in d = 3, the integral over the wavenumber diverges logarithmically at small scales. 
Therefore, the Lenard-Balescu equation is marginally valid provided that a small-scale cut-off is introduced at the 
Landau length Xmin — X^ (sec the previous argument). In the thermal bath approach, the diffusion tensor is given 
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by Eq. (11-39) with^ 



2^3(x) = — -T^4 / ' _ .... rffc- (20) 



2(27r)2 n + S(a;)fc2,]^ + C(a;)2fc 
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This integral converges at large scales (fc — > 0) due to Debye shielding. For small velocities |v| — > 0, the diffusion 
coefficient is given by (11-45) with^ 



Vf^" = 2?3(0) = ^^—k% I ^ dk. (22) 

' ^ 2(27r)i « ""io (k^ + klf 

This corresponds to the Dcbyc-Hiickcl approximation. In the dominant approximation In A 3> 1, the Lenard-Balescu 
diffusion coefficient (|20p and Debye-Hiickel diffusion coefficient reduce to the Landau diffusion coefficient 
with InA = ln(AD/Ami„) in which the Debye length appears naturally. This shows that the Landau equation with a 
large-scale cut-off at the Debye length provides a very good approximation of the Lenard-Balescu equation in d = 3. 
This implies that dynamical screening is ncgligible^'^ in d = 3. Indeed, in the dominant approximation, wc can neglect 
the velocity effects encapsulated in the dielectric function e(k, k • v) and use the Debye-Hiickel potential accounting 
for static screening. 

For a Coulombian plasma in d = 3 (marginal case), the impact theory and the wave theory are two approximate 
theories that arc complementary to each other. Collective interactions between charges are not included in the impact 
theory implying that the integrand in the Boltzmann equation is valid only for impact parameters sufficiently smaller 
than the Debye length (A <C Xd)- On the other hand, the curvature of orbits at small impact parameters is not 
included in the wave theory implying that the integrand in the Lenard-Balescu equation is valid only for wavelengths 
sufficiently larger than the Landau length (A ^ Al). For A ^ 1. the range of validity of these two theories greatly 
overlaps in the region Al <C A <C Ad corresponding to the domain of validity of the Landau equation. A possibility 
to unify these approaches and determine precisely the numerical factor in the Coulomb logarithm was first recognized 
by HulDbard [l^l and subsequently developed by Kiliara and Aono [s^ and Gould and DeWitt [33| among others. 
Roughly speaking, the idea is to sum the Boltzmann (0 < A < Ad) and Lenard-Balescu (Al < A < +oo) equations 
and subtract the Landau equation (Al < A < Xd)- This leads to a fully convergent kinetic equation. It has the 
form of the Landau equation in which the Coulomb logarithm is exactly given by In A — 27 -I- In 2 where 7 = 0.5772 
is Euler's constant. If we neglect the constant term as compared to In A (dominant approximation), this procedure 
justifies using the Landau equation with a small-scale cut-off at the Landau length and a large-scale cut-off at the 
Debye length. 

Conclusion: The evolution of the system as a whole is described by the Lenard-Balescu equation (II-3) with a cut-off 
at the Landau length (this cut-off is not ad hoc but is justified by the correct treatment of strong collisions). The 
relaxation of a test particle in a thermal bath is described by the Fokker-Planck equation (11-36) with the diffusion 
tensor given by Eqs. (11-39) and (PO)) . In the dominant approximation In A ^ 1, the Lenard-Balescu equation can be 
replaced by the Landau equation (II-5) with a cut-off at the Debye length (this cut-off is not ad hoc but is justified 
by the correct treatment of collective effects). In that case, the evolution of the system as a whole is described by the 
Landau equation (II-5) or (II-9) with K3 = /m^) ^(Ad/Al). The relaxation of a test particle in a thermal bath 
in described by the Fokker-Planck equation (11-36) with the diffusion tensor given by Eqs. (11-47) and (11-48) where 
I?3 given by Eq. (jlSp . Using the results of Sec. IIIIl we find that the relaxation time is given by 

Tfl^v'^ A 

~ / 3 \ 1 — I^D, (23) 

ne4 1n(!2£!|.) InA 



3 



The value of the integral is 



Jo 



= dk = \nk+ - In 



B(x)2 + C(x)2 



B{x) 
2\C{x)\ 



■K B(x) 
arctan 

2 \C{x)\ 
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The value of the integral is 



fk 

Jo 



fc3 



■ dk : 



ln(l + A^) 



A2 



1-I-A2 



(19) 



(21) 



{k^ + klf 

^'^ This is true for the Maxwell distribution. This is not true anymore for distributions for which the large velocity population plays a 
crucial role as compared to the bulk of the distribution. 
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where to ^ ojp is the dynamical time. A short historic of the early development of the kinetic theory of 3D 
Coulombian plasmas with many references is given in Paper II. 



B. 2D Coulombian plasmas 



If we neglect strong collisions and collective effects, the evolution of the system is described by the Landau equation 
(II-5). In the thermal bath approach, the diffusion tensor is given by Eq. (11-45) with 

This integral converges at small scales implying that strong collisions arc negligible. However, it diverges rapidly 
(linearly) at large scales implying that collective effects (Debye shielding) are important. If we introduce a large-scale 
ko at the Debye length, we obtain 

^2 = -^vlkol. (25) 



2(2^)t " ^A' 



where A = nk^ ^ 1 represents the number of electrons in the Debye disk. 

Collective effects are taken into account in the Lenard-Balescu equation (II-3). This equation is rigorously valid for 
a Coulombian potential in c? = 2 (see Sec. IIIIB[) . The large-scale divergence that occurs in the Landau equation is 
regularized by Debye shielding. In the thermal bath approach, the diffusion tensor is given by Eq. (11-39) with 

1 „3 /• + 00 .2 

I?2(x) = T—kh 5 dk. (26) 

2(27r)i n ^ 7o [k^^ + B{x)klf + Cix^k^ 

This function is perfectly well-defined. For small velocities |v| — )• 0, the diffusion coefficient is given by Eq. (11-45) 
with 

^ ^'"^ - r = = i^-'- 

This corresponds to the Debye-Hiickel approximation. For |v| -> -l-oo, we can use the asymptotic results given in Sec. 
IV. D of Paper I. The diffusion coefBcient D^{v) is given by Eq. (II-53-d) with 2?2 replaced by I?2(0) and the diffusion 
coefficient D\\{v) is given by Eq. (II-53-c) with I?2 replaced by 2?2(1) where 

I vt A f^°° k^ 1 o 1.45804 , , 

2?2 1 = - — kj) / 5 dk = -vf„kD 7 = 1.45804I?2. 28 

2(2^)1 n [P + S(l)fcy' + C(l)2fc4 2(2^)1" A 

We note that the Dcbye-Hiickel approximation is not correct at large velocities (compare Eqs. (|27|) and ((28)) ). 
Therefore, dynamical screening is important for 2D plasmas, contrary to 3D plasmas in the dominant approximation. 
This makes the dimension d — 2 particularly interesting. On the other hand, the Lenard-Balescu diffusion coefficient 
(PS)) and the Debye-Hiickel diffusion coefficient (P7)) are different from the Landau diffusion coefficient (PS)) . Therefore, 
the Lenard-Balescu equation cannot be rigorously approximated by the Landau equation with a large-scale cut-off at 
the Debye length. However, comparing Eqs. ([25)) . ()27p and (|25)) . we see that the discrepancy is not dramatic since 
the velocity dependence of the diffusion coefficient is the same (at least asymptotically) and the value of the prefactor 
is only slightly different (7r/4 = 0.785398 and 1.45804 instead of 1). Therefore, the Landau equation may be used as 
a rough approximation of the Lenard-Balescu equation by "adapting" the value of the large scale cut-off. 

Conclusion: The evolution of the system as a whole is described by the Lenard-Balescu equation (II-3). The 
relaxation of a test particle in a thermal bath in described by the Fokker-Planck equation (11-36) with the diffusion 
tensor given by Eqs. (11-39) and (|25)) . The Landau equation (II-5) with a cut-off at the Debye length is not rigorously 
equivalent to the Lenard-Balescu equation in d = 2 but the discrepancy is not dramatic. If we use this approximate 
theory, the evolution of the system as a whole is given by the Landau equation (II-5) or (II-9) with K2 = 2iTe^ /m?kD, 
or = -K^e^ /2m^kD in the Debye-Hiickel approximation. The evolution of a test particle in a thermal bath is 

described by the Fokker-Planck equation (11-36) with the diffusion tensor given by Eqs. (11-51) and (11-52) where I'2 
given by Eq. (|25)) . or by Eq. (|27| in the Debye-Hiickel approximation. Using the results of Sec. IIIIl we find that the 
relaxation time is given by 

- ^^Wx ~ ^^D: (29) 
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where ~ is the dynamical time. 

A kinetic theory of 2D Coulombian plasmas has been developed by Benedetti et al. [s^ for the cut-off Coulombian 
potential and by Chavanis [s^ for the true Coulombian potential. However, collective effects are treated in an 
approximate manner in Appendix C of that paper. The present treatment is more satisfactory. 



C. ID Coulombian plasmas 

If we neglect strong collisions and collective effects, the evolution of the system is described by the Landau equation 
(II-5). In the thermal bath approach, the diffusion coefficient is given by Eqs. (11-46) and (11-54). This leads to 

This integral converges at small scales implying that strong collisions arc negligible. However, it diverges rapidly 
(quadratically) at large scales implying that collective effects (Debye shielding) are important. If we introduce a 
large-scale cut-off fc„iin ~ kjj at the Debye length, we obtain 

- 2(2,)./'(ff„.)V.A '-'°-''. 

where A = nkjj^ ^ 1 represents the number of electrons in the Debye segment. 

Collective effects are described by the Lenard-Balescu equation (IT3). This equation is rigorously valid for a 
Coulombian potential in d = 1 (see Sec. IIII Bp . The large-scale divergence that occurs in the Landau equation is 
regularized by Debye shielding. In the thermal bath approach, the diffusion coefficient is given by Eqs. (11-43) and 
(11-44) leading to 

1 3 r+oo IF 

T I. ^ (32) 



B{^/^v]kl 



D 



CiJf^v) k% 



The integral is perfectly well-defined and can be calculated analytically. Using the identity 



+00 



*i(y) = (^2 _^^)2 _^ I " I ^ ^ arctan(2/) = ^arcot(y), (33) 



we obtain 



D{v) = ^ . arcot I | . (34) 

27r(/3m)^A|z;| I /^|^;|e-5/3™-"' 

Since arcot(.T) ^ 1/x for x — > +00 and arcot(a;) — > tt for x — >■ —00, we have the asymptotic behaviors 

^(^^ " 2(2^)V2(^™)3/2A' ^(") 2{PmYK\v[ ^^^^ 

We note that, for v = Q. the exact diffusion coefficient given by Eq. ([35l -a). coincides with the expression (|3T|) obtained 
from the Landau equation by introducing a large-scale cut-off at the Debye length. However, this is essentially 
coincidental. Indeed, the exact diffusion coefficient for v = Q may be written as 



1 v: 



3 r+00 
'A 



m = 77^ — k% jr^^dk. (36) 



(27r)i/2 n (k^ + kl)^ 

This corresponds to the Debye-Hiickel approximation. It turns out that the value of this integral is the same as 
the value of the integral ([5(1)) with a large-scale cut-off at kn. However, this is not generally true. Collective effects 
are not equivalent to introducing a large-scale cut-off at ko in the Landau diffusion coefficient, even for v = (see, 
e.g.. Sec. IIVB[) . On the other hand, the diffusion coefficients ([51]) and are totally different at large velocities 
\v\ — )■ +CX1 [compare Eqs. ([51]) and ([551 -b)]. Indeed, the exact diffusion coefficient ([M]) decays algebraically as 
while the diffusion coefficient (|3ip based on the Landau equation with a large-scale cut-off at the Debye length decays 
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exponentially rapidly as e"'^™" 1"^ . This is a crucial difFcrencc. This clearly shows that, in d = 1, the Landau equation 
with a cut-off at the Debye length gives wrong results (except for low velocities). Collective effects must therefore be 
properly taken into account by using the Lenard-Balescu equation. 

Conclusion: The evolution of the system as a whole is described by the Lenard-Balescu equation (ITS). Actually, 
this equation reduces to df /dt = in d = 1 (see Eq. (ILll)). Therefore, the collision term vanishes at the order 1/A. 
If it does not vanish at the order l/A^, the relaxation time of the system as a whole is given by 

^whole ^ ^rn _ ^2^^ > ^^^^ (37) 

The relaxation of a test particle in a thermal bath is given by the Fokker-Planck equation (n-42) with the diffusion 
coefficient Using the results of Sec. IIIIl we find that the relaxation time of a test particle in a thermal bath is 

given by 

tT" ^ ^ Md, (38) 

where ^ ujp^ is the dynamical time. We note that the relaxation time is independent on the density. 

As discussed in Paper II, in d = 1 the relaxation time of a test particle in a thermal bath is different from the 
relaxation time of the system as a whole. Since the distribution of the field particles does not change on a timescale 
of the order Ai^i, we can consider the relaxation of a test particle in an out-of-equilibrium bath (see Sec. IV. C of 
Paper II). This relaxation process is described by the Fokker-Planck equation (11-63) with a diffusion coefficient given 
by Eq. (11-58) and (11-61). Using the identity ([55)1 . it can be written as 

= '?iy\| arcot (-^^V ^ dv) . (39) 

7rm|/'(w)| V Afim J-oo u-v J 

For the isothermal distribution, we recover Eq. (|34|) . For the watcrbag distribution, using Eq. (11-64), we get 

2' 



D[v) 



2m 



1 - 



(40) 



for -~v,n < V < and D{v) = otherwise. When collective effects are neglected, the diffusion coefficient is given 
by Eq. (11-65). However, it involves a divergent integral as in Eq. (|30| for the isothermal distribution, so this 
approximation is not reliable. 

A one component ID plasma model with a neutralizing background^^ has been studied numerically by Dawson 
pol l4l| . This is one of the first models in physics where numerical simulations have been made to test non-equilibrium 
properties of particles in interaction. The kinetic theory of this model, taking collective effects into account, has been 
developed by Eldridge and Feix [i^, [l^^^^ what concerns the evolution of the system as a whole, they found that 
the Lenard-Balescu collision term vanishes in ID so that the relaxation time is strictly larger that Ato (numerical 
simulations show that it scales like A^t^ which corresponds to the next order in the expansion of the BBGKY 
hierarchy). On the other hand, considering the relaxation of a test particle in a thermal bath, they found that the 
evolution of the distribution function is governed by the Fokker-Planck equation (11-42) with the diffusion coefficient 
In that case, the relaxation time scales like i^*'* ~ Ato- Recently, the Dawson model has been re- investigated 
by Sano and Kitahara [1^. The present results complete their theoretical analysis. 



V. SELF-GRAVITATING SYSTEMS 



For self-gravitating systems, the Fourier transform of the potential of interaction is 



{2nru{k) = m = ^, (41) 



The equilibrium properties of a two-components ID plasma where the electron and ions move freely has been studied by Eldridge and 
Feix [37II . This corresponds to the Lenard-Prager model [ssl. |39| . 
We independently obtained the same results in Sec. 2.8.3 of |44| . 
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and the Fourier transform of the Debye-Hiickel potential is 

{27rrunH{k)^~j^, ^nH{k)^j^. (42) 

fC J J 

In physical space, we have u{r) ~ —G/r and UDnir) = —Gcos{kjr)/r in d = 3, u{r) = Glnr and uonir) ~ 
^GYo{k,jr) in d = 2, u{x) = G\x\ and udh{x) ~ (G/fcj) sin(fcj|x|) in d = 1. 



A. 3D self-gravitating systems 

If we neglect strong collisions, collective effects, and make a local approximation, the evolution of the system is 
described by the Vlasov-Landau equation (11-16). In the thermal bath approach, the diffusion tensor is given by Eq. 
(11-45) with 

V, ^ ^^^fc^ / ^. (43) 
2(2^)5 n ^ k 

This integral diverges logarithmically at small and large scales. This means that both strong collisions and spatial 
inhomogcneity must be taken into account. However, their effect is weak since the divergence is only logarithmic (we 
arc in the marginal case of Sec. IIII Cp . If we introduce a small-scale cut-off k^ax = = Afc,/ at the Landau length 
at which binary collisions become strong and a large-scale cut-off k„iin = kj at the Jeans length which represents the 
typical size of the system, we get 

"-i^iW""'''^' '''' 

where A = nk^^ ^ N ^ 1 represents the total number of stars in the cluster. 

Strong collisions are taken into account in the Boltzmann equation. This corresponds to the two-body encounters (or 
impact) theory. In this equation, the integral over the impact parameter does not diverge at small scales. Therefore, 
the correct treatment of strong collisions regularizes the logarithmic divergence that appears at small scales in the 
Landau equation. However, for a Newtonian potential in d = 3, the integral over the impact parameter diverges 
logarithmically at large scales. Therefore, the Boltzmann equation is marginally valid for a Newtonian potential in 
d = 3 provided that a large-scale cut-off is introduced at the Jeans length Xmax = Aj (see the next argument). Since 
weak collisions dominate over strong collisions, we can expand the Boltzmann equation for small deflexions (or directly 
use the Fokker-Planck equation). In the dominant approximation InA 3> 1, we obtain the Landau equation (11-16) 
with the diffusion coefficient (1441) with InA = hi{Xmax/XL) in which the Landau length Al appears naturally (see 
Appendix [D| . 

Spatial inhomogeneity and collective effects are taken into account in the Lenard-Balescu equation written with 
angle-action variables [l6l - l20| . This corresponds to the wave theory. If we neglect collective effects, we get the Landau 
equation written with angle-action variables. In these equations, the integral over the wavenumber does not diverge 
at large scales since the finite extent of the system is taken into account. Therefore, the correct treatment of spatial 
inhomogeneity regularizes the logarithmic divergence that appears at large scales in the Landau equation. However, for 
a Newtonian potential in d = 3, the integral over the wavenumber diverges logarithmically at small scales. Therefore 
the Lenard-Balescu and Landau equations written with angle-action variables are marginally valid for a Newtonian 
interaction in d = 3 provided that a small-scale cut-off is introduced at the Landau length Xmin = A^ (see the previous 
argument). If we neglect collective effects, make a local approximation, and introduce a large scale cut-off at the Jeans 
length, we get the Vlasov-Landau equation (11-16) with the diffusion coefficient (|44p with In A = ln(Aj/Amira). In the 
dominant approximation In A ^ 1, this equation does not depend on the precise value of the large-scale cut-off since it 
appears in a logarithmic factor. This suggests that the local approximation is reasonable and that the Vlasov-Landau 
equation (II-3) provides a good approximation of the true dynamics of a stellar system in d = 3. 

For a stellar system in d = 3 (marginal case) , the impact theory and the wave theory are two approximate theories 
that are complementary to each other. Collective interactions between stars are not included in the impact theory 
implying that the integrand in the Boltzmann equation is valid only for impact parameters sufficiently smaller than 
the Jeans length (A ^ Aj). On the other hand, the curvature of orbits at small impact parameters is not included in 
the wave theory implying that the integrand in the Lenard-Balescu equation with angle-action variables is valid only 
for wavelengths sufficiently larger than the Landau length (A ^ Al). For A ^ 1, the range of validity of these two 
theories greatly overlaps in the region A^ <C A ^ Aj corresponding to the domain of validity of the Vlasov-Landau 
equation. Combining these approaches, we can motivate (but not rigorously justify) using the Vlasov-Landau equation 
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with a small-scale cut-off at the Landau length and a large-scale cut-ofF at the Jeans scale. A better kinetic equation 
is the Lenard-Balescu equation written with angle-action variables with a small-scale cut-off at the Landau length. 

Conclusion: The evolution of the system as a whole is described by Lenard-Balescu equation (taking collective 
effects into account) or by the Landau equation (neglecting collective effects) written with angle-action variables with 
a small-scale cut-off at the Landau length |16H20| . The evolution of a test particle in a thermal bath is described by 
the Fokker-Planck equation written with angle-action variables with a small-scale cut-ofF at the Landau length. If 
we neglect collective effects, make a local approximation, and introduce a large scale cut-off at the Jeans length, the 
evolution of the system as a whole is described by the Vlasov-Landau equation (11-16) or (11-17) with ~ 27rmG^ In A. 
The relaxation of a test particle in a thermal bath is described by the Fokker-Planck equation (11-56) with the diffusion 
tensor given by Eqs. (11-47) and (11-48) where is given by Eq. (|33])^^. Using the results of Sec. we find that 
the relaxation time is given by 



tR 



nm?G'^ In 



/2„3/2|33/2 



N 



(45) 



where tjj ^ is the dynamical time. Collective effects tend to reduce the value of the relaxation time as explained 
in Appendix [CI A short historic of the early development of the kinetic theory of 3D stellar systems with many 
references is given in [20j . 



B. 2D self-gravitating systems 

If we neglect strong collisions, collective effects, and make a local approximation, the evolution of the system is 
described by the Vlasov-Landau equation (11-16). In the thermal bath approach, the diffusion tensor is given by Eq. 
(11-45) with 

^2 - —m — kj / TT- (46) 



2(27r)i n ' Jo P 

This integral converges at small scales implying that strong collisions are negligible. However, it diverges rapidly 
(linearly) at large scales implying that spatial inhomogeneity effects are important. If we introduce a large-scale 
cut-off fcmm = kj at the Jeans length, we obtain 

where A = nfcy^ ~ ^ 1 represents the number of particles in the system. 

Spatial inhomogeneity and collective effects are taken into account in the Lenard-Balescu equation written with 
angle-action variables |16l - [2(j |. This equation is rigorously valid for a Newtonian potential in c? = 2 (see Sec. IIIIBp . 
If we neglect collective effects, we get the Landau equation written with angle-action variables. These equations 
converge at large scales since they take into account the finite extent of the system. If we make a local approximation, 
and introduce a large scale cut-off at the Jeans length, we get the Vlasov-Landau equation (11-16). Contrary to the 
3D case, the local approximation is not very accurate in 2D because the divergence at large scales in the diffusion 
coefficient (|46p is linear instead of logarithmic. Therefore, the value of the diffusion coefficient is strongly dependent 
on the large-scale cut-off. However, the Vlasov-Landau equation (11-16) may provide a reasonable approximation of 
the true dynamics of self-gravitating systems in c? = 2 by properly "adapting" the value of the large-scale cut-off. 

Conclusion: The evolution of the system as a whole is described by Lenard-Balescu equation (taking into account 
collective effects) or by the Landau equation (ignoring collective effects) written with angle-action variables [T6l - [20j . 
The evolution of a test particle in a thermal bath is described by the Fokker-Planck equation written with angle-action 
variables. If we neglect collective effects, make a local approximation, and introduce a large-scale cut-off at the Jeans 
length, the evolution of the system as a whole is described by the Vlasov-Landau equation (11-16) or (11-17) with 
K2 = 2TiG^m/kj. In the thermal bath approach, we get the Fokker-Planck equation (11-56) with the diffusion tensor 



For self-gravitating systems, a pure isothermal distribution /(r, v) = Ae~'^'"(" /2+*(r)) ig^ds to a density profile with an infinite mass. 
Therefore, a pure isothermal distribution function cannot hold in the whole cluster. We can, however, apply the Fokker-Planck theory 
in the core of the system which is approximately isothermal. 
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given by Eqs. (11-51) and (11-52) where given by Eq. (1471) -'-^. Using the results of Sec. IIIIl we find that the 
relaxation time is given by 

^« ^ G3/2„l72^3/2 - ^^D, (48) 

where to ^ <^7^ is the dynamical time. Collective effects tend to reduce the value of the relaxation time as explained 
in Appendix [Cl 

In a recent paper, Marcos [49| developed a kinetic theory of self-gravitating systems in two dimensions. He showed 
through extensive numerical simulations that a local approximation may be implemented in 2D gravity provided that 
a large-scale cut-off is introduced (it is a fraction of the system's size). He also derived the diffusion and friction 
coefficients from an extension of the Chandrasekhar binary collision theory. This leads to the expressions (II-B6- 
a) and (II-B6-b) written in terms of the Rosenbluth potentials (II-B7). He then simplified these expressions for a 
Maxwellian distribution leading to his equations (7a) and (7b). We previously 0,llll obtained the simpler expressions 
(11-35), (H-51) and (11-52) from the Landau approach. As explained in Appendix [Pl the two approaches (which make 
a local approximation and neglect collective effects) are equivalent. The collisional relaxation of 2D self-gravitating 
systems has also been studied by Teles et al. [50| using a simplified dynamics. 



C. ID self-gravitating systems 

If we neglect strong collisions, collective effects, and make a local approximation, the evolution of the system is 
described by the Vlasov-Landau equation (11-16). In the thermal bath approach, the diffusion coefficient is given by 
Eqs. (11-46) and (11-54) leading to 

1 ii3 , o /■ + °° rlh 



This integral converges at small scales implying that strong collisions are negligible. However, it diverges rapidly 
(quadratically) at large scales implying that spatial inhomogeneity effects are important. If we introduce a large-scale 
cut-off k,nin = kj at the Jeans length, we obtain 

D(v) = — g-^P^y^ f50) 

^ ' 2(27r)i/2(/3„j)3/2A^ ' ^^"^ 

where A = nkj^ ^ N ^ 1 represents the number of particles in the system. 

Spatial inhomogeneity and collective effects are taken into account in the Lenard-Balescu equation written with 
angle-action variables |16l - [20| . This equation is rigorously valid for a Newtonian potential in d = 1 (See Sec. IIII Bp . If 
we ignore collective effects, we get the Landau equation written with angle-action variables. These equations converge 
at large scales since they take into account the finite extent of the system. If we neglect collective efi^ects, make a local 
approximation, and introduce a large scale cut-off at the Jeans length, we get the Vlasov-Landau equation (11-16). 
However, the Vlasov-Landau equation is expected to display important discrepancies with the true dynamics for at 
least two reasons, (i) First of all, it is clear that we cannot make a local approximation to describe the evolution of 
the system as a whole. Indeed, in d = 1 the collision term in the Vlasov-Landau equation (11-16) vanishes so this 
equation reduces to the Vlasov equation. Therefore, it implies that the relaxation time scales as N'^to (like for ID 
plasmas) while numerical simulations [53| of ID self-gravitating systems show that the relaxation time is of order 
Nt£). Actually, when spatial inhomogeneity is properly accounted for, additional resonances appear in the Lenard- 
Balescu equation (or in the Landau equation) written with angle-action variables, and the collision term is non-zero 
(see discussion in [l3|)- This explains why the relaxation time scales as Ntjj instead of N'^to- (ii) On the other 
hand, in the thermal bath approach, the Vlasov-Landau equation (11-16) leads to the Fokker-Planck equation (11-57) 



In d = 2, the velocity dispersion of a self-gravitating system in a steady state is exactly given by (-u^) = 2v'^ = GM/2 (this result can 
be derived from the virial theorem [26l . l4tl ) . Statistical equilibrium states exist at a unique temperature ksT = GMm/A. On the other 
hand, the density profile of a distribution of field particles at statistical equilibrium, and the gravitational potential that it generates, are 
known analytically (see, e.g., [47l. |48|| and references therein). They are given by p(r) = po/(l-|-7rpor^/M)^ and '3?(r) = ''^('Jf^ "'"''^) 
where po is the central density of the cluster which parameterizes the series of equilibria. These expressions may be substituted in the 
local diffusion coefficient II47I I and in the advection term of the Fokker-Planck equation (11-56). 
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with the diffusion coefficient ([50]) "'^^. Contrary to the 3D case, the local approximation is not valid in ID because the 
divergence at large scales in the diffusion coefficient is quadratic instead of logarithmic. Therefore, the value of 
the diffusion coefficient is strongly dependent on the large-scale cut-off. Furthermore, we should not give too much 
credit on the precise expression ([50]) of the diffusion coefficient. We have seen in the case of ID plasmas that collective 
effects change the velocity dependence of the diffusion coefficient (for large velocities). This is probably true also for 
ID self-gravitating systems. For all these reasons, the local approximation is not a good approximation in c? = 1. 

Conclusion: The evolution of the system as a whole is described by Lenard-Balescu equation (taking collective 
effects into account) or by the Landau equation (neglecting collective effects) written with angle-action variables 
p^| - [20| . The evolution of a test particle in a thermal bath is described by the Fokker-Planck equation written with 
angle-action variables. Using the results of Sec. Illli we find that the relaxation time is given by 

ifl ~ ^ ~ NtD, (51) 
Gm 

where to ^ ijJq^ is the dynamical time. We note that the relaxation time is independent on the density. Collective 
effects tend to reduce the value of the relaxation time as explained in Appendix [CJ The local approximation cannot 
be employed for ID self-gravitating systems. 

The collisional relaxation of ID self-gravitating systems has been studied by Miller [5l|, Valageas [Hj, Joyce and 



Worrakitpoonpon [53| and Sano [5 



VI. THE HMF MODEL 



The HMF model 55| consists in N particles moving on a circle and interacting via a potential u{9 — 9') = -^[1 — 
ecos(0 — 9')] with e = +1 in the attractive case and e = — 1 in the repulsive case. The Fourier transform of the 
potential is m„ = jjg^{2Sn,o — £Sn,±i) and its normalized Fourier transform is T]n = ■^l3{—2dn,o + ef^n.ii)- The repulsive 
HMF model does not display any phase transition. Its statistical equilibrium states are always spatially homogeneous 
(for any energy and temperature). The attractive HMF model displays a second order phase transition. Its statistical 
equilibrium states are spatially homogeneous for E > Ec = 3/4 (i.e. T > Tc = 1/2) and spatiallyinhomogeneous for 
E < Ec (i.e. T < Tc); when E < Ec (i.e. T < Tc), the spatially homogeneous phase is unstable [j]. Here, we restrict 
ourselves to the stable homogeneous phase (the kinetic theory of the spatially inhomogcneous HMF model can be 
treated with angle-action variables as in [l6|-|20|). 

If we neglect collective effects, the relaxation of a test particle in a thermal bath is described by the Fokker-Planck 
equation (11-42) with a diffusion coefficient given by Eqs. (11-46) and (11-54) leading to 



D(v) = e-3P^ . (52) 



If collective effects are taken into account, the diffusion coefficient is given by Eqs. (11-43) and (11-44) leading to 

^^7^ • (53) 

he -13 + V2/33/2i;e-5^^" /V 2" ^y^^y \ _^ p3„2g-/3t,^ 



N 



It has the asymptotic behaviors 



TV (2e-;3)2' ' ' N2 \27r 

The exact expression (|53p of the diffusion coefficient reduces to the approximate expression (|52l) at high temperatures 
T +00 or for large velocities \v\ — > +00. For low velocities v ^ 0, the diffusion coefficient ([5^ a) can be obtained 



In d = 1, the density profile of a distribution of field particles at statistical equilibrium, and the gravitational potential that it generates, 
are known analytically (see, e.g., |47|| and references therein). They are given by p{x) = {GM'^m./AkBT)cosh-^{GMmx/2kBT) and 
$(3;) = {2kBT /m)\a[cosh{GMmx/2kgT)]. These expressions may be substituted in the local diffusion coefficient I I50II and in the 
advection term of the Fokker-Planck equation (11-57). 
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from Eq. ((53)) by replacing the "dressed" potential 
M„/|e(n, 0)|. In physical space, it reads 



dressed 



udh( 



1 

N 



1 



1-^6/3 

+1, the diffusion coefficient £'(0) diverges as (2 



|e(n, by the Debye-Hiickel potential u„ 



') 



DH 



(55) 



In the attractive case e = +1, the diffusion coefficient £'(0) diverges as (2-/3) ^ at the critical point 
Introducing the dynamical time to ~ ^n/vm, we find that the relaxation time t^*'' ^ v^/D scales like 



bath 



e?Nt, 



(56) 



where we have estimated the diffusion coefficient at i; = (i.e. in the Debye-Hiickel approximation). This formula 
shows that the relaxation time depends on the temperature. At high temperatures, we can neglect collective effects 
and we obtain i^''*'' - AT'^Nto- On the other hand, if we take collective effects into account, we find that the 
relaxation time tends to zero at the critical point Tc = 1/2 in the attractive case e = +1. Therefore, collective 
effect reduce the relaxation time close to the critical point. The same conclusion has been reached for self-gravitating 
systems (see Appendix [C]). However, this prediction could be tested numerically more easily with the HMF model 
since the potential of interaction is restricted to one Fourier mode and the system is finite. 

The evolution of the system as a whole is given by the Lenard-Balescu equation (II-3) which reduces to df /dt ~ 
for spatially homogeneous distributions in d = 1 (see Eq. (II-ll)). This implies that the relaxation time of the system 
as a whole is 



uhole 



> Nt, 



(57) 



If the system remains always spatially homogeneous, we expect that the relaxation time scales as N'^tj), like for ID 
plasmas (see Sec. IIV C[) . due to the absence of resonances at the order 1/A^. If the system remains always spatially 
inhomogeneous, we expect that the relaxation time scales as Nto, like for ID stellar systems (see Sec. IV C[) . due to 
additional resonances at the order 1 /N brought by spatial inhomogeneity. Finally, if the system is initially spatially 
homogeneous but, due to the development of correlations (finite N effects), becomes Vlasov unstable and experiences a 
dynamical phase transition from a homogeneous phase to an inhomogeneous phase [H^ , we expect that the relaxation 
time is intermediate between the two previous scalings Ntu and N^tp. This argument is consistent with numerical 
results showing that the relaxation time scales as N^-'^tp in that case 
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Since the distribution of the field particles does not change on a timescalc of the order Nt o in the homogeneous 
case, we can consider the relaxation of a test particle in an out-of-equilibrium bath (see Sec. IV. C of Paper II). This 
relaxation process is described by the Fokker-Planck equation (11-63) with a diffusion coefficient given by Eqs. (11-58) 
and (11-61). When collective effects are neglected, we obtain 



When collective effects are taken into account, we get 
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(59) 



If f{v) decreases sufficiently rapidly, we find that the exact diffusion coefficient ([59]) behaves as D{v) ~ N^^tt"^ f{v) 
for \v\ ^ +00 as when collective effects are neglected. For the isothermal distribution (thermal bath), we recover Eq. 
(j53p . For polytropic (Tsallis) distributions, the diffusion coefficient D{v) has been computed numerically in [ssj . For 
the waterbag distribution, using Eq. (11-64), we get 



D{v) 



1 IT 



(60) 



if —Vm < V < Vm and D{v) ~ otherwise. As noted in [S^, for the attractive HMF model, the diffusion coefficient 
diverges when the velocity v of the test particle is equal to the pulsation uj = {v^ 

-1/2)1/2 Qf^jjg 

wave arising from the 

slightly perturbed distribution of the field particles. This divergence occurs for any stable distribution {vm > l/V^ i.e. 
E > Ec) because, for the waterbag distribution, the modes are purely oscillatory. For a single humped distribution 
with a maximum at v — vq, the diffusion coefficient ((5^ diverges only at the critical point {E = Ec) for v — vq. 
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Indeed, the critical point E = Ec is such that e(l, vq) = [s^. For the isothermal distribution, this implies that D{0) 
diverges when T — >■ in agreement with Eq. (|54l-a). 

The kinetic theory of the HMF model was first considered by Inagaki [g^I- He tried to adapt the Lenard-Balescu 
equation to that model but, unfortunately, his paper contains mistakes that led him to incorrect conclusions. The 
relaxation of a test particle in a thermal bath was considered by Bouchet [6l] | who derived the diffusion coefficient 
([53|) from the study of the stochastic process of equilibrium ffuctuations. In later works by Bouchet and Dauxois [HH 
and by Chavanis 0, , it was realized that the kinetic theory of the HMF model could be performed by adapting 
the results of plasma physics to this specific system. Indeed, the previous results can be viewed as particular cases of 
the general results presented in Paper II. 



VII. CONCLUSION 



In this series of papers, we have completed the literature on the kinetic theory of systems with long-range inter- 
actions. A general formalism has been developed in Papers I and II and it has been applied to specific potentials of 
interaction in the present paper. We have considered pure power-law potentials of interaction and we have explained 
how the divergences that occur in the Landau equation may be cured by taking into account strong collisions (Boltz- 
mann) or collective effects (Lenard-Balescu) depending on the value of the index 7. We have also treated the case of 
plasmas and stellar systems in d dimensions and the case of the HMF model with attractive or repulsive interactions. 
The relation to previous works has been discussed in detail. An interest of our presentation is to provide a "unified" 
picture of the subject. 



Appendix A: The different kinetic equations 

The standard kinetic equations (Boltzmann, Fokker- Planck, Vlasov, Landau, Lenard-Balescu) may be derived from 
the Klimontovich equation by using a quasilinear approximation [l^ as discussed in Paper I. They may equivalently 
be derived from the Liouville equation by neglecting three-body correlations in the BBGKY hierarchy [§, [ll|, [l^l ■ In 
this Appendix, writing the BBGKY hierarchy in a symbolic form, we show the connection between these different 
equations and discuss their domains of validity without going into technical details. 

For spatially homogeneous systems, the first two equations of the BBGKY hierarchy may be written symbolically 

as 

% = C[9], ^ + {Co+C')g+C[f,g]+nh]=S[f], (Al) 

where / is the one-body distribution function, g the two-body correlation function, and h the three-body correlation 
function. In the first equation, the collision term C[g] describes the effect of two-body correlations on the evolution of 
the distribution function. In the second equation, C = Cq + C is a two-body Liouvillian operator where Co describes 
the free motion of the particles and C describes the exact two-body interaction. The term C [/, g] describes collective 
effects and the term T[h] describes three-body correlations. Finally, S[f] is a source term depending on the one-body 
distribution function. 

For systems with short-range interactions, the potential decreases at large distances as with 7 > d (i.e. a > 2). 
In general, there is also a repulsion at short distances like in the Lennard- Jones potential or in the hard sphere 
potential. Therefore, the potential has a finite range Iq. In the dilute limit uIq ^ 1 (this corresponds to <C 0; 
can neglect collective effects and three-body collisions (C = T = 0). Therefore, the evolution of the system is driven 
by weak and strong collisions. The coupled equations (jAip with C = T ~ lead to the Boltzmann equation fl^ . 

For systems with long-range interactions, the potential decreases at large distances as r~'^ with 7 < d (i.e. a < 2). 
In the the weak coupling approximation 9 ^ 1 or A = nXj^ 3> 1 (this corresponds to Xjj ^ I), at the order 1/A, 
we can neglect three-body correlations (T = 0) and strong collisions (£' = 0). Therefore, we can replace the exact 
two-body interaction of the particles by their free motion {Cq + — Cq) which amounts to making a linear trajectory 
approximation. The evolution of the system is driven by weak collisions and collective effects. The coupled equations 
(|Aip with C' ~T = lead to the Lenard-Balescu equation If, in addition, we neglect collective effects (C = 0), 
we get the Landau equation^^. 



We must be careful, however, that for potentials that are singular at r = 0, strong collisions may be important at small scales (in that 
case, the expansion of the BBGKY hierarchy in powers of g = 1/A is not uniformly convergent). This may invalidate the Lenard-Balescu 
and Landau equations. This difficulty has been discussed in Sec. Illll for purely power-law potentials. 
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The Landau equation (£' = 7" = C = 0) is intermediate between the Boltzmann and the Lenard-Balescu equation. 
It describes the effect of weak colUsions but ignores strong colhsions and collective effects. It can be obtained from 
the Lenard-Balescu equation (£' = T = 0) by neglecting collective effects (C = 0) or from the Boltzmann equation 
(C = T = 0) by considering the limit of small deflexions [C = 0). 

Actually, these kinetic equations describe the effect of collisions at different scales (the scale A may be interpreted 
as the impact parameter). For A ^ Al (small impact parameters), the collisions are strong and we must solve the two- 
body problem exactly and take into account the bending of the trajectories of the particles. For A I (intermediate 
impact parameters), the collisions are weak and we can make a weak coupling, or straight line, approximation. For 
A A^i (large impact parameters), we must take collective effects into account. The Boltzmann equation is valid 
for A ^ Xd- It describes strong collisions (A ^ A^) and weak collisions (A ~ I). The Lenard-Balescu equation is 
valid for A ^ A/,. It describes weak collisions (A ~ /) and collective effects (A ~ A^i). The Landau equation is valid 
for Al ^ A <C \d- It describes weak collisions. When we go beyond the domains of validity of these equations, 
divergences occur, and appropriate cut-offs must be introduced. 

For spatially inhomogeneous systems (see [l^l for more details), the first two equations of the BBGKY hierarchy 
may be written symbolically as 

|[ + (Vo + V,„./.[/])/ = ^ (£0 + £' + [/])5 + C[/, 5] + T[h] = S[f]. (A2) 

In the first equation, V — Vq + ^m.f. is the Vlasov operator taking into account the free motion Vq of the particles 
and the advection by the mean field Vm./.- In the second equation, the term Cm.f. in the two-body Liouvillian takes 
into account the effect of the mean field in the two-body problem. 

For systems with long-range interactions, in the weak coupling approximation 5 ^ 1 or A = nXj 3> 1 (this 
corresponds to Aj » /), at the order 1/A, we can neglect three-body correlations (T = 0) and strong collisions 
(£' = 0). Therefore, we can replace the exact two-body dynamics of the particles by their mean field motion 
(£0 + C + Cm.f. — C-o + Cm.f.)- The evolution of the system is driven by weak collisions, collective effects, and 
spatial inhomogeneity. The coupled equations (jA2p with £' = T = lead to the Lenard-Balescu equation written 
with angle-action variables |16l - [20| ^^. If, in addition, we neglect collective effects (C = 0), we get the Landau equation 
written with angle-action variables. Finally, if we neglect collisions (which is rigorously valid for A — > -|-oo) we get 
the Vlasov equation. 

Appendix B: The modified Landau equation 

For pure power-law potentials with a > ac (i.e. 7 > 7c or d > 3 for the Coulombian or Newtonian potential), the 
diffusion coefficient ([3|) is divergent at small scales (k -l-oo) while it is convergent at large scales (fc — >■ 0). In that 
case, it is necessary to take strong collisions into account, while collective effects (and spatial inhomogeneity effects 
for attractive interactions) may be neglected (see Sec. IIII Ap . The dynamical evolution of the system is then described 
by the Boltzmann equation. However, for long-range potentials (a < 2, i.e. 7 < d), weak collisions dominate over 
strong collisions when A 1. We can therefore expand the Boltzmann equation for weak deflexions |Av| <C I. This 
is equivalent to starting from the Fokker-Planck equation (11-25) and computing the first and second moments of 
the velocity increments by a binary encounter theory. However, in the calculation of (Avi) and (AviAvj), we must 
use the exact trajectory of the particles (and take into account the strong deflexions for collisions with small impact 
parameters) in order to have convergent expressions. This leads to the modifled Landau equation, written below in 
d = 3 for a > ttc = (i.e. 7 > 7c = 1): 

^ = A _d_ r , wH,,-w,Wj 

dt m\m) dvij ^ u'(3-")/(i+«) 

where ^4 is a coefficient (depending on a) related to the integration over the impact parameter. For a > ac, the 
modified Landau equation (|Bip does not present any divergence (see Sec. IIII Ap . For a < ac, the modified Landau 
equation ()Bip presents an algebraic divergence at large scales and it is not valid anymore. In that case, it must 
be replaced by the Lenard-Balescu equation (II-3) as explained in Sec. IIIIBI For a = ac, Eq. (|Bip reduces to 
the ordinary Landau equation (II-9). This equation presents a logarithmic divergence at large scales which can be 
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^-^)/(v,0/(v',0, 



(Bl) 



We can use angle-action variables because for A ^ 1 the coUisional evolution is a slow process which allows us to make an adiabatic 
approximation. 
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cured by introducing a large-scale cut-off at the Debye length Xmax = (properly justified) as explained in Sec. 
nil CI On the other hand, Eq. (|Bip does not diverge at small scales since the effect of strong collisions is taken into 
account. In that case, the Landau length appears naturally in the diffusion coefficient = {27Te'^ /ni^)\n{Xmax/ ^l) 
(see the discussion in Appendix ID|) . We emphasize that, for a > ac, the modified Landau equation (jBip differs from 
the ordinary Landau equation (II-9) since the power of the relative velocity w appearing in the denominator of the 
collision term is (3 — a)/(l + a) instead of 3 (in d = 3 dimensions). From Eq. (|B1|) we find that the relaxation time 
scales as 

3-ct 

tR- ""I ^A^tn, (B2) 

in agreement with Ec^. ([6l). 

The modified Landau equation (jBip was introduced by Potapenko et al. [fi^l- The Coulombian case corresponds 
to 7 = 1 and the case of Maxwell molecules corresponds to 7 = 4. The potential is soft for 1 < 7 < 4 and hard for 
7 > 4. Some applications will be discussed elsewhere. 



Appendix C: Reduction of the relaxation time due to collective effects for attractive interactions 

For attractive power-law potentials of interaction with a < Oc, the evolution of the system is described by the 
Lenard-Balescu equation written with angle-action variables |16l420| . This equation takes spatial inhomogeneity and 
collective effects into account. Unfortunately, it is extremely complicated. If we make a local approximation, we 
obtain the Vlasov-Lenard-Balescu equation (11-15). However, this equation presents a strong divergence at the Jeans 
scale. Indeed, in the thermal bath approximation, the diffusion tensor is given by Eq. (11-39) with 



2(27r)'^-5 A Jo [k^-^~B{x)]' + C{xy 



T^<ii^) = TTTTTTTTrr^ / . . „. ^, .,2 . dn. (Cl) 



If we make the Debye-Hiickel approximation, or consider small velocities |v| — > 0, the diffusion tensor is given by Eq. 
(11-45) with 

nDH 1 "m^^ f~^°° '^'^ I,. ,bath A 



2(2^)<^-^ A 7o {K^^-^-l? ' " /+-^_^^ 

Clearly, the integral diverges at k — kj (i.e. k = 1). This is of course related to the Jeans instability for a s patially 
homogeneous system. This divergence does not occur when spatial inhomogeneity is properly accounted for |16l42C| . 
This divergence suggests, at a heuristic level, that collective effects (which account for anti-shielding) tend to increase 
the diffusion coefficient, hence to decrease the relaxation time. This reduction should be particularly strong for a 
system close to instability due to an enhancement of fluctuations [13, ■ These arguments are further developed in 
[2I Hi] and in the Appendix E of [1^. 

These results show that it is not possible to make a local approximation and simultaneously take collective effects 
into account. The usual procedure is to ignore collective effects, make a local approximation, and introduce a large- 
scale cut-off at the Jeans length. This is the procedure that is usually followed in stellar dynamics jl3l4l5l. A more 
rigorous approach is to use the Lenard-Balescu or Landau equations written with angle-action variables |16l420| . 



Appendix D: Connection between the Landau and the Chandrasekhar kinetic theories 

In this Appendix, we discuss the relation between the Landau kinetic theory and the Chandrasekhar kinetic theory. 
These kinetic theories were developed independently and their connection was not fully realized by early workers on 
the subject (see footnote 5 in Paper II). Even in modern textbooks of astrophysics [l^ll^: the kinetic theory of stellar 
systems is presented with the approach of Chandrasekhar and the Landau equation is not mentioned. 

Landau [66[ developed a kinetic theory of 3D Coulombian plasmas by starting from the Boltzmann equation and 
using a weak deflexion approximation. In his calculations, he approximated the trajectories of the particles by straight 
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lines even for collisions with small impact paramcters"'^^. This leads to the kinetic equation (II-9) with 

j^Landau = ^ In ('^'j . (Dl) 
fn \ \min J 

This factor presents a logarithmic divergence at small and large scales that Landau regularized heuristically by 
introducing a small-scale cut-off at the Landau length {\min = ^l) and a large-scale cut-off at the Debye length 
i^max = Ad). 

Chandrasekhar [631 developed a kinetic theory of 3D stellar systems by starting from the Fokker-Planck equation 
(II-25)-(II-26) and calculating the first and second moments of the velocity increments resulting from a succession 
of binary collisions. In the calculation of {Avi) and (AviAvj), he used the exact trajectory of the stars (i.e. he 
solved the two-body problem exactly) and took into account the strong deflexions due to collisions with small impact 
parameters'^. In the dominant approximation InA^ ^ 1, his approach completed by Roscnbluth et al. [ssj leads 
to the expressions (II-B6) of the diffusion tensor and friction force expressed in terms of the Roscnbluth potentials 
(II-B7) with 

j^Chandra ^ ^^^^2 ^"g^^m ^ ^2) 

This factor presents a logarithmic divergence at large scales that can be regularized heuristically by introducing a 
cut-off at the Jeans length^°. However, contrary to the Landau approach, this factor does not present a logarithmic 
divergence at small scales since the effect of strong collisions is taken into account explicitly in the Chandrasekhar 
approach^'. As a result, the gravitational Landau length Xl = Gm/v^ appears naturally in the calculations of 
Chandrasekhar. Of course, these results can be translated to the plasma case by replacing G by (e/m)^ and A/ by 
Xd- 

The Landau and the Chandrasekhar kinetic theories make the same assumption: binary encounters and an expansion 
in powers of the momentum transfer. They actually differ in the order in which these are introduced. Landau starts 
from the Boltzmann equation and considers a weak deflexion approximation while Chandrasekhar directly starts from 
the Fokker-Planck equation but calculates the coefhcients of diffusion and friction with the binary-collision picture. At 
that level, their theories are equivalent since the Fokker-Planck equation can precisely be obtained from the Boltzmann 
equation in the limit of weak deflexions. The crucial difference is that Landau makes a weak coupling assumption 
and ignores strong collisions while Chandrasekhar takes them into account. 

Except for this important difference, the calculations of Appendix B of Paper II show that the kinetic equation (II-9) 
derived by Landau [6^ is equivalent to the kinetic equation (II-B8) derived by Chandrasekhar [63] and Roscnbluth 
et al. [68j] although they appear under a different form. This equivalence is not always apparent in the astrophysical 
literature (see the discussion in [2^ ) . In the astrophysical literature, the diffusion and friction coefficients are usually 
derived from the Chandrasekhar binary encounter theory leading directly to Eqs. (II-B6)-(II-B7) while they can be 
obtained equivalently from the Landau equation (II-Bl) leading to Eqs. (II-B3)-(II-B4) which can be transformed 
into Eqs. (II-B6)-(II-B7) as we have shown in Appendix B of Paper II. 

We stress that the equivalence between the Chandrasekhar approach and the Landau approach is true for any 
dimension d > 1 and for any long-range potential of interaction (their limitation has been discussed in this paper) . 
Therefore, the kinetic theory developed by Gabrielli et al. [2l| and Marcos [i^ based on the Chandrasekhar approach 
is equivalent to the kinetic theory that we have presented in [s^ based on the Landau approach. Indeed, essentially 



This is equivalent to starting from the Fokker-Planck equation (11-25)- (11-26) and calculating the first and second moments of the 
velocity increments resulting from a succession of binary collisions by using a straight line approximation. This is also equivalent to 
starting from the BBGKY hierarchy and neglecting three-body collisions, collective effects, and strong collisions (T = C = £' = 0) as 
discussed in Appendix [XI 

This is equivalent to expanding the Boltzmann equation for weak deflexions while taking the effect of strong collisions (i.e. the bending 
of the trajectories) into account. This is also equivalent to starting from the BBGKY hierarchy and neglecting three-body collisions and 
collective effects {T = C = 0), but retaining strong collisions {C 7^ 0), as discussed in Appendices lAl and iBl 

Actually, Chandrasekhar introduced a large-scale cut-off at the inter-particle distance I but this choice was later on criticized by many 
authors. It is now acknowledged that the proper large-scale cut-off in a stellar system is the Jeans length Aj which is the gravitational 
analogue of the Debye length Xjy in plasma physics. 

The approach of Chandrasekhar takes into account strong collisions with an impact parameter smaller than the Landau length Aj^ 
which yield a deflection at an angle larger than 90°. As we have seen, this can suppress the divergence at small scales. However, the 
Chandrasekhar approach docs not take into account the possibility of forming binary stars which correspond to bound states with strong 
correlation between particles (see Sec. V of Paper II). In other words, Chandrasekhar only considers hyperbolic trajectories and not 
elliptical ones in the two-body problem. The effect of binary stars in the kinetic theory requires a special treatment [69j . 
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the same assumptions are made: spatial homogeneity (or local approximation), neglect of collective effects, and weak 
coupling approximation. The Lenard-Balescu approach that we have presented in [7| is more general since it can take 
collective effects into account. Finally, the generalized Landau approach presented in [ITl . 1201 can take into account 
spatial inhomogeneity and the generalized Lenard-Balescu approach presented in [l8l.ll9l| can take into account spatial 
inhomogeneity and collective effects. 



Appendix E: The choice of the large-scale cut-off 

The kinetic theories developed by Landau (66| in plasma physics and by Chandrasekhar (67| in stellar dynamics, 
modeling the collisional process by a succession of independent weak binary encounters, lead to a diffusion coefficient 
and a friction force that diverge logarithmically at large impact parameters. These authors heuristically circumvented 
this problem by introducing a large scale cut-off. 

In the case of plasmas, Persico [tO], Landau [iBi, Bohm and AUer [7l|, Prigogine and Balescu [tH, and Cohen et 
al. [7^ argued that the logarithmic divergence should be cut-off at the Dcbye length . This leads to a Coulombian 
factor 



where A = nA^, is the number of electrons in the Debye sphere. The kinetic theory of Lenard and Balescu proved 
that the Debye length indeed represents the relevant large-scale cut-off to introduce in the Landau equation (this is 
exact in the dominant approximation). 

In stellar dynamics, the choice of the large-scale cut-off Xmax created some debate in the early literature. Jeans 
[23| . Spitzer [tI], Chandrasekhar and von Neumann [t^, and Prigogine and Balescu [t^] argued that the logarithmic 
divergence has to be cut-off at the interparticle distance I ^ n~^/^ (see Sec. Ill Cp . This leads to a Coulombian factor 

\n(^\^\n(l-\=l\nK, (E2) 



where A = nAj ~ iV is the number of stars in the Jeans sphere. However, Cohen et al. [t^ argued that the divergence 
has to be cut-off at the Jeans scale Aj (of the order of the system's size) which is the gravitational analogue of the 
Debye length. This leads to a Coulombian factor 

The corresponding relaxation time is lower than the original Chandrasekhar relaxation time by a factor 3/2. Hcnon 
[7^, Prigogine and Sevcrne [77| argued that if the system were spatially homogeneous and infinite, the relevant 
large-scale cut-off would be the mean free path A (sec Sec. IIIIDp . This leads to a Coulombian factor 



In ( = In A ) ^ 2 In A. (E4) 



The corresponding relaxation time is lower than the original Chandrasekhar relaxation time by a factor 3. Finally, 
Severne and Haggerty [t^ showed that when the spatial inhomogeneity and the finite extent of the system are properly 
taken into account, there is no large-scale divergence anymore. This suggests that the Jeans scale is the most relevant 
large-scale cut-off to consider. 
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